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I. INTRODUCTION 

The properties of point defects in metal oxides and 
their superlattices have been an area of recent interest, 
and have been investigated using increasingly precise and 
sophisticated experimental methods. [1, 2] Optical exper- 
iments can determine differences between sharp photo- 
luminescence peaks with meV precision; [3] this can be 
used to determine the different types of defects, but the 
chemical identity of a point defect or defect complexes 
remains difficult to obtain. Electronic structure methods 
have been used to obtain the formation energies and the 
location of defect levels within the gap to identify the 
source the photoluminescence peaks. [4, 5] However, first 
principles calculation of defects in complex systems suf- 
fer from several limitations. [6] The long-range interaction 
between defects in neighboring cells, an artifact of peri- 
odic boundary conditions (PBC), can affect relaxation 
around a defect. [7] Also, correcting the spurious inter- 
action between charged defects [8, 9] and the introduced 
compensating background charge is an area of ongoing 
interest. Furthermore, there is a well-known underesti- 
mation of band gaps when calculations using the local 
spin density approximation [10] (LSDA) of density func- 
tional theory (DFT) are performed; while it is possible 
to apply a post-hoc correction to the gap and get good 
results for uniform solids, [11] it is less obvious how com- 
parable to experiment the result of this practice is when 
applied to the location of the defect levels. More ad- 
vanced semilocal density functionals perform better but 
still generally underestimate the band gap. [12] 

Some of these limitations may be overcome via post- 
DFT or hybrid DFT calculations, which often give band 
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gaps closer to experiment. [12] Band gap accuracy gen- 
erally implies efficacy in modeling defect structures; for 
example, a recent comparison of semilocal and hybrid 
density functional theory (DFT) calculations showed 
that there is a strong correlation between the calcu- 
lated valence band width (VBW) and defect formation 
energies. [13] Accurate VBW values obtained with the 
Heyd- Scuseria- Ernzerhof screened hybrid functional [14] 
(HSE) lead to point defect formation energies and energy 
levels in close agreement with experiment for elemental 
and binary nonmetals. [15-17] These findings are encour- 
aging because (unlike semilocal functionals) screened hy- 
brid functionals require no post-DFT correction of the 
band gap and provide excellent lattice parameters in ad- 
dition to other bulk properties. [18-21] In studies of cu- 
bic SrTiOs (STO) [22-25] and other pervoskites,[23, 26- 
28] HSE has been shown to perform well for band gaps, 
spin-phonon effects, and numerous structural effects, in- 
cluding strain, bulk moduli, octahedral angles tilts and 
rotations amongst others. Recently, a regular and hy- 
brid density functional study on many perovskites [29] 
has highlighted the relevance of anharmonic corrections 
to lower the octahedral rotation and tilt angles. The an- 
harmonic correction was 0.15° for all compounds except 
STO for which the correction was 0.8°. For methods 
containing 25% HF mixing, namely PBEO and PBEsolO, 
the uncorrected angles in SrTiOs were already close to 
experiment; applying the anharmonic correction lead to 
too small angles. HSE belongs to the same family but has 
a screening parameter uj = 0.11, so it can be thought of 
as an interpolation between PBE (u =oo) and the global 
hybrid PBEO (uo =0). We expect the anharmonic effects 
correction would lead to HSE octahedral rotations over- 
corrected to too small angles compared to experiment. 

We here present the results of our HSE calculations 
for the formation energies and band structures of neu- 
tral defects in SrTiOs. This work is motivated by HSE's 
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agreement with experiment in the calculation of many of 
the electronic, structural and elastic properties in cubic 
SrTiOa[22, 23] and other metal oxides. [26] It is worth not- 
ing that the direct and indirect band gaps as well as VBW 
are in excellent agreement with experiment (see table I) 
compared to the results from LSDA and B3PW [30], 
which were previously used to study point defects in 
SrTi03. We focus on the various neutral vacancies in 
SrTi03, such as the oxygen (Vq), strontium (Vg r ) and ti- 
tanium vacancies (V^), as well as on the effect of doping 
STO with La (La 5r ). 

II. COMPUTATIONAL METHODS 

All calculations included in this manuscript were per- 
formed using the development version of the GAUSSIAN 
suite of programs, [31] with periodic boundary conditions 
(PBC)[32] used throughout. The functionals applied in 
this work include the generalized gradient approxima- 
tion of Perdew, Burke and Ernzerhof[33, 34] (PBE) and 
HSE.[14] The Def2-[35] series of Gaussian basis sets were 
optimized following our procedure, described in Ref. 22. 
We use the notation SZVP to differentiate these opti- 
mized PBC basis sets from the molecular Def2-SZVP ba- 
sis sets. Strontium has the inner-shell electrons replaced 
with smallcore pseudopotentials, while for oxygen and 
titanium atoms we utilize all electron basis sets. SZVP 
basis offers a good compromise between computational 
efficiency and high accuracy for electronic and structural 
properties of bulk SrTi0 3 . [22] The use of SZVP basis set 
with HSE screened hybrid functional (HSE/SZVP) im- 
poses limitations on the size of the supercell that can be 
efficiently computed, so a STO supercell of 2x2x2 (40 
atoms) was used with a dense /c-point mesh of 6x6x6 
which included the T point. Calculations with larger su- 
percells of 2x3x3 (90 atoms) with the same density of 
/c-points were performed in order to discuss the impor- 
tance of defect self-interactions and the effect of varying 
the defect concentration on the electronic properties of 
STO. The pruned integration grid for DFT employed was 
(99,590), which corresponds to the Gaussian option "ul- 
trafine". The remaining numerical settings in GAUSSIAN 
were left at the default values, e.g., geometry optimiza- 
tion threshold was set to 450 x 10 -6 hartrees/bohr, SCF 
convergence was set to "tight". Unless otherwise noted, 
crystal structures for chemical potential calculations on 
SrO, TiO, Ti203, and Ti20 were downloaded as CIF files 
from the ICSD.[36] All structures are fully relaxed (opti- 
mized) unless otherwise noted. 



III. RESULTS AND DISCUSSION 

A. Vacancy Formation Energy Calculations 

The calculations of neutral defect formation energies 
are based upon the formalism of Zhang and Northrup.[41] 



TABLE I. Comparison of properties of bulk STO relevant 
to defect formation from our work, previous computational 
studies, and experiment. Our HSE/SZVP band gaps and va- 
lence band widths (VBW) are closer to experiment than the 
semilocal LSDA and the global hybrid B3PW results. 





LSDA a 


B3PW a 


HSE b 


Exp. 


Direct gap (eV) 


2.36 


3.96 


3.59 


3.75 c 


Indirect gap (eV) 


2.04 


3.63 


3.20 


3.25 c 


VBW (R-»R)(eV) 


4.77 


6.47 


5.0 


5.0 C 


ao(A) 


3.86 


3.90 


3.89 


3.89 d , 3.90 e 


B(GPa) 


214 


177 


204 


179 



a Ref. 37 
b Ref. 22 
c Ref. 38 
d Ref. 39. 
e Ref. 40. 



These values are generated using the equation 
E f = E T - [^(perfect) 

TlSrf^Sr ^Tif^Ti ^of^o] 

where E T and £r (perfect) are the calculated total en- 
ergies of the supercells containing the point defect and 
the perfect bulk host materials, respectively. The num- 
ber of each element removed from the perfect supercell 
is represented by n x , while ji x corresponds to the atomic 
chemical potentials in an SrTiOs crystal. Given the as- 
sumption that SrTiOa is always stable, the chemical po- 
tentials of the these elements can vary with the following 
restriction: 

^Sr + I^Ti + 3/io = ^SrTi0 3 (bulk) ( 2 ) 

Atomic chemical potentials vary according to the sam- 
ple composition and cannot be determined exactly. How- 
ever, they can be varied to cover the whole phase diagram 
of SrTi0 3 , by splitting into SrO, TiO, Ti 2 3 , and Ti 2 
bulk phases. Hence, the calculated formation energies 
for the neutral point defects vary according to equilib- 
rium positions; an example of this is O-rich versus O- 
poor points on the phase diagram. The calculated en- 
thalpies of formation in idealized materials (non-relaxed 
structures) for phases containing Sr, Ti and O are sum- 
marized in Table II, and are compared to previous LSDA 
calculations [42] and experiment. [43] 

As a general trend, the formation enthalpies computed 
with HSE are close to the results from semilocal func- 
tionals LSDA and PBE (from this work), although the 
HSE values are slightly higher. The only exception is 
SrO, where both HSE and PBE tend to overestimate the 
formation enthalpies to the same extent, exceeding the 
LSDA values. 

The formation energies of vacancy defects in STO as 
function of its composition are plotted in Figure 1 with 
the points A to G [44] based on the phase diagram in 
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TABLE II. Calculated enthalpies of formation in eV/atom 
for idealized materials for phases containing Sr, Ti and O 
compared to previous LSDA calculations and experiment. 





HSE 


PBE 


LSDA a 


Exp. b 


Ti0 2 


-3.92 


-3.75 


-3.76 


-3.24 


Ti 2 3 


-3.74 


-3.60 


-3.63 


-3.15 


TiO 


-2.99 


-2.95 


-3.04 


-2.81 


SrO 


-4.00 


-3.90 


-3.36 


-3.07 



a Ref. 42 
b Ref. 43 





25 r 




20 - 


'(eV) 


15 - 


D) 




CD 
C 




CD 


10 - 


C 

o 




Format 


5 - 
- 



vo 

VTi 
VSr 



A B C D E F G 

FIG. 1. (Color online) Defect formation energies of isolated 
neutral vacancies in STO at each equilibrium point based 
upon the phase diagram detailed in the text. 

Ref. 42. The plot shows that HSE's formation energies 
show a similar dependence on the chemical potentials as 
was observed for the LSDA results of Ref. 42. However, 
both V^i and Vg° r curves are shifted to higher formation 
energies, while the Vq curve is slightly shifted downwards 
compared to the LSDA spectrum. At point A, Vg° r has 
the lowest formation energy followed by and Vq. At 
point B, Vg° r remains the most stable but Vq becomes 
more stable than V^ { contrary to the LSDA predictions. 
Under low partial pressures of oxygen (O-poor conditions 
corresponding to points C to G) , Vq vacancies form more 
easily and dominate the spectrum while the formation 
energies of V^ { and Vg° r keep increasing. 

A quantitative comparison with previously published 
calculations using the same cell size formation ener- 
gies [30, 42, 45] available at point B (oxygen rich con- 
ditions) reveal several issues of interest. 

• For Vq, the semilocal functional PW91[45] yields 
a formation energy of 8.56 eV, the global hybrid 
B3PW [30] predicted a higher formation energy of 
8.74 eV, while LSDA [42] gave a lower value at 
7.95 eV. Our HSE value is the lowest among these, 
yielding a formation energy of 7.43 eV. This de- 
crease in the Vq formation energy is one of the fac- 
tors leading to a strong competition with V® { (see 
below). 



• For Vy' V the semilocal functional LSDA [42] yields 
a formation energy of 5.7 eV making it far more 
stable than Vq. In contrast, with HSE we predict 
that Ef = 7.86 eV meaning that V® { becomes less 
stable than Vq. 

• Vg° r remains the most stable defect under these con- 
ditions, with HSE providing a value of 2.81 eV com- 
pared to the 1.7 eV from LSDA. 

Overall, the quantitative formation energies differences 
between our HSE results and previous LSDA results are 
substantial. These might originate from the enhanced ac- 
curacy we gain with HSE in the calculated valence band 
width (VBW) and band gap of bulk SrTi0 3 (see table I) 
identified in Ref. 13 to reflect an enhanced precision in 
the defect formation energies. 

It is also worth mentioning size effects. The forma- 
tion energy of Vq at point B in a larger supercell of 90 
atoms, giving a 1.8% Vq concentration, is 7.66 eV, which 
is 0.2 eV higher than the value obtained with the smaller 
supercell. 




FIG. 2. (Color online). Ideal, crystalline SrTiOs, as observed 
from the 100 face, where the oxygen anions are in red, the 
strontium cations are in green and the titanium cations are 
in blue. Magenta, yellowish and cyan circles highlight the 
locations of the isolated Vq, Vs° r and V^i defects, respectively. 



B. Isolated Vacancies: Relaxation Effects 

Isolated, neutral defects have been introduced into the 
crystal structure of cubic STO by removing one atom of 
either O, Sr or Ti, respectively, as depicted in Figure 2. 
The structure was then fully relaxed using HSE/SZVP. 
Figure 3 illustrates the major displacements in each of 
the defect structures, namely Vq, Vg° r and V^- v relative to 
the idealized crystal. Distance decreases are depicted by 
black arrows, while increasing lengths are denoted by yel- 
low arrows. In each case, the magnitude of displacement 
is implied by the size of the arrows: larger /thicker arrows 
indicate greater deviation from the defect-free structure. 

A 4% Vq vacancy was created by removing one oxygen 
(Figure 3-a) from a 2x2x2 cell, inducing small, asym- 
metric bond length changes throughout the cell and an 




FIG. 3. (Color online). Relaxed SrTiOs with a single isolated vacancy. Each structure is in the 100 orientation and the location 
of the vacancy is noted. The distortions induced by Vq, Vs° r and V^i are highlighted (a), (b) and (c), respectively. Oxygen 
anions are red, Sr cations are green, and Ti cations are blue. 



overall decrease in volume relative to the crystal. At the 
vacancy, the remaining O-Ti bonds all contract but to 
different extents. Near the vacancy, oxygens move to- 
ward nearby strontium, inducing O-Sr bond contractions 
of up to 5.66 pm. Further from the vacancy, in the re- 
gions of the cell behind the defect (and the plane of the 
figure), strontium appears to migrate inward, toward the 
vacancy, causing the more distant O-Sr bonds to lengthen 
by as much as 6.05 pm. 

The Vq structure contracts, and undergoes a tetrago- 
nal distortion where the lattice parameters are a— 3.879 
and 6=c=3.890 A, compared to the 3.902 A in the crystal. 
High concentration of oxygen vacancies were identified to 
be responsible of a similar tetragonal distortion reported 
experimentally [3] in STO doped with La. Worth not- 
ing that with HSE/SZVP, we are able to capture this 
small tetragonal distortion corresponding to a c/a ratio 
of 1.0028. This distortion was also observed in previous 
LDA calculations [46] for supercells of the same size, but 
the c/a ratio was not reported. 

As seen from Figure 3-b, the atomic relaxation around 
the V$ Y induce staggered bond lengths (vertically) and 
a general contraction of O-Sr bonds, with the large de- 
crease being 1.55 pm, smaller than those observed for O- 
Sr in the presence of an oxygen vacancy. The alternating 
distances cause the linear Ti-O-Ti angle to decrease to 
173.9°. Again, the defect structure contracts, with all 
lattice parameters decreasing: a=6=c=3.880A . 

The V^i structure depicted in Figure 3-c shows the 
axial oxygens experiencing a 9.2 pm decrease in length, 
while the medial oxygens move away from the vacancy by 
0.6 pm. The Sr-0 distances shrink over a range of 1.51 to 
5.40 pm. While also asymmetric, these distortions have 
less of an effect on the cell., with the equilibrium lattice 
parameters being a=c=3.892 and 6=3.894 The a/b ratio 
here is very small so the structure remains effectively 
cubic following full relaxation. 

In summary, the loss of metallic species results in 



smaller volumes with the loss of the larger metal, Vg° r , 
producing greater contraction. Removing the oxygen 
produces larger deviations from cubic symmetry in ad- 
dition to overall shrinkage of the unit cell. 

C. Isolated Vacancies: Band Structures 

The electronic band structure and the projected den- 
sities of states (PDOS) of the ideal and nonideal/defect 
STO supercells (doped with O, Sr and Ti neutral vacan- 
cies) are shown in figure 4. To allow a clear comparison 
between the different systems, the VBM was set as ref- 
erence for all systems while the Fermi energy is depicted 
with a solid line. The PDOS are plotted alongside with 
their corresponding band structure and rescaled in the 
same way. 

For the ideal 2x2x2 STO supercell the Brillouin zone 
folds, causing a direct band gap of 3.19 eV. The CBM 
in the ideal solid is triply degenerate at T and composed 
of a heavy electron band (he), a light electron band (le) 
and a spin-orbit band (so). The band degeneracy is lifted 
as we move far from the V point towards high symmetry 
directions X, M and R, giving rise to different values 
of the electron effective masses ra^, mie, mso in each 
direction. The electronic effective mass for each band i 
can be computed by the curves of energy versus k being 
fitted to a parabola [47] ; k is taken from the T point up 
to 0.5% along r X, T -> M and T R paths using 
the formula: 

1=2^ (3) 
m H dk 2 

We estimate the T — >• X effective masses, without 
spin-orbit correction, vnF h ^ x and mf^ x to be 7.3m e 
and 0.5m e respectively, where m e is the free electron 
mass. This agrees well with previous nonrelativistic 
calculations [23, 47, 48] using LSDA and HSE, which 
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K-Path PDOS (arb . imits) K_Path PDOS (arb. onite) 



FIG. 4. (Color online). Band structures and PDOS calculated with HSE/SZVP for the 2x2x2 SrTiOs supercell. The top 
figures represent bulk STO, and V§ while the bottom row contains Vs° r and VSi- The Fermi energy Ef is indicated by a solid 
black line. The red bands indicate the occupied defect bands while the unoccupied defect bands are shown in blue. 



predicted he masses ranging from 6.1 to 7.3m e and 
le masses of 0.4m e . The remaining effective masses 
are ra^ M =l.lra e , mf^ M =0.8m ei m^ R =0.9m e and 
mf^ R =0.7m e . 

Introducing into STO, the defect band averaged in 
the Brillouin zone is located at 0.13 eV above the valence 
band maximum (VBM) and 3.06 eV below the conduc- 
tion band maximum (CBM) (see Table III) . The defect 
band has a VB character and is triply degenerate at T 



TABLE III. Location of the various defect bands aver- 
aged in the Brillouin zone with respect to the conduction 
band minimum ( A6db-cbm ) and valence band maximum 
(A6db-vbm). Data are for the 2x2x2 SrTiOs supercell com- 
puted with HSE/SZVP. Defect formation energies (in eV) at 
point B (oxidation condition) are reported as well. 





AeDB-CBM 


AeDB-VBM 


E/ 




3.06 


0.15 


7.86 




3.27 


0.07 


2.81 


vS 


0.44, 0.40 a 


2.82, 2.90 a 


7.43 



a Experimental estimations from Refs. 3 and 49 



with a single band occupied by two electrons, while the 
upper two bands are empty (see figure 4). The crys- 
tal field resulting from the atomic relaxation around 
causes the triply degenerate CBM band to split in the 
following way: the heavy electron band (he) is followed 
by the light electron band (le) at 26 meV, while the spin- 
orbit band (so) is located 14 meV higher. 

For F s ° r , a triply degenerate defect band appears at 
0.07 eV above the VBM and 3.27 eV below the CBM 
at the T point. From figure 4, two bands are occupied 
by four electrons while the upper band is empty. The 
degeneracy of the defect band is lifted at the X point, 
where the two occupied band remain degenerate while 
the unoccupied band is 600 meV higher in energy. The 
he and le bands in the CBM remain degenerate while the 
so band is located at 1.5 meV higher in energy mainly 
because the overall deformation V$ r introduces into the 
lattice is small. 

Introducing Vq into STO causes a defect band (DB) 
populated with two electrons to appear in the gap just 
below the bulk-like conduction band maximum labelled 
here as (DCBM), indicating a donor band. The Fermi 
energy is shifted from the top of the valence band to 
the maximum of the defect band at the M point fol- 
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lowed by the empty conduction band. The CBM con- 
serves the bulk character as it remains empty and triply 
degenerate (clearly shown in Figure 4) corresponding 
to the Ti t2# states followed by a doubly degenerate e g 
band. The above CBM conservation of degeneracy con- 
firms that the band appearing underneath is a fully oc- 
cupied non-degenerate defect band. The isosurface of the 
highest occupied state at the T special k-point (Figure 5) 
shows that the electronic charge density is localized on 
Vq occupying the Ti dangling bonds. 

We calculated the position of the defect state by av- 
eraging over the Brillouin zone. We found that the DB 
is located in average at 0.44 eV below the CBM for the 
40 atoms supercell (see table III) and 0.42 eV for the 90 
atoms supercell. This could be compared the recent ex- 
perimental measurements, which place the position of the 
defect level 0.4 eV below the CBM, causing the blue light 
photoluminescence of STO at room temperature. [3, 49] 
In the present calculations, self-defect interaction or the 
so-called size effects were minimized by keeping the de- 
fects neutral, and conserving the same high density of 
fc-point sampling of the Brillouin zone. [50] Nevertheless, 
simulations using larger supercells and finite size scal- 
ing [7] are still needed for a more accurate comparison. 
As of today's computational resources, a full relaxation 
of larger supercells using HSE/SZVP and the computa- 
tional settings used here are computationally very expen- 
sive. This defect band seems to be shallower than the one 
computed with B3PW, which was located 0.79 eV below 
the CBM. [30] We attribute this to the B3PW functional 
overestimating the band gaps (indirect gap of 3.6 eV com- 
pared to 3.2 eV from experiment). Our HSE-predicted 
defect level for Vq is deeper than that calculated with 
LDA; in those calculations the defect band can often not 
be distinguished from the CBM (resonant band), [46] al- 
though in some cases it lies as little as 0.08 eV below the 
CBM. [42] 

The atomic relaxation around Vq leads to a tetragonal 
distortion of the supercell where c/a= 1.0028 (see Sec- 
tion IIIB), corresponding to a —0.3% compressive strain 
along the X axis. This amount of strain is comparable to 
the values applied experimentally in La doped STO [24] 
and identified as causing a substantial increase of the 
electron mobilities (by a factor of 3.3) due to the appear- 
ance of light electron effective mass along the strain and 
the transport direction. 

In our calculation, the effective mass m\^ x /m e along 
the direction parallel to the compressive strain (fc" = 
7r/a(100)) drops from the bulk value of 7.3 to 0.5 in de- 
fective SrTi03(se table IV), indicating that the electron 
mobility increases substantially along the x-axis. How- 
ever, m^ x /m e increases along the k 1 - = 7r/a(010) = 
7r/a(001) directions (perpendicular to the compressive 
strain) to 11.0 indicating a reduced electron mobility in 
these directions. Along r —> M, the effective massesre- 
main unchanged in the fc" = 7r/a(110) = 7r/a(101) con- 
taining the strain axis. 



D. STO doped with La 

The substitution of one Sr atom by a La atom in our 
2x2x2 STO supercell lead to a dopant concentration of 
12.5%, which is low enough to be compared to experi- 
mental data, [3] where dopant concentrations as high as 
15% has been used. The resulting Sr .875Lao.i25Ti03 
compound relaxes to a cubic structure with lattice pa- 
rameters a=6=c=3.893 A. This agrees well with experi- 
ments [3] which demonstrated that La doping conserves 
the cubic symmetry, and has a negligible effect on the 
c-axis lattice parameter. 

Figure 6 shows the band structure of 
Sro.875Lao.i25Ti03 for spin up and spin down elec- 
trons. For the spin up electrons, the VBM is almost 
triply degenerate with a very small splitting of 7 meV; 
this is also the case of the lowest CB, a band which is 
populated with one extra electron and shows a small 
band broadening (13 meV). The energy difference 
between the last VB and the first populated CB is 3.16 
eV. The next set of conduction bands, normally located 
at 200 meV above the CBM in the bulk STO, is split 
into one band a 172 meV above the CBM followed 
by a doubly degenerate band 100 meV higher. The 
visualization of isosurface of the highest occupied state 
at the T special k-point show that the charge density is 
not localized on the defect site, but rather corresponds 
to the Ti t2 9 states. The electronic density of states 
(not shown) also confirms that CB remains dominated 
by Ti 3d states as it is does in the STO bulk phase; La 
Ad starts to contribute to the CBM only at about 2.3 eV 
above the Fermi energy. 

For the spin down electrons (bottom of figure 6) the 
band structure conserves most of the bulk band charac- 
teristics with triply degenerate CBM and VBM, but the 
band gap is larger than the bulk value, measuring 3.55 
eV. The heavy electron effective masses relative to the 



TABLE IV. Non-relativistic HSE/SZVP calculated heavy 
electron effective masses (in units of the free-electron mass 
m e ) at the conduction maximum along different high- 
symmetry directions in ideal (m^ BM ) and oxygen vacancy 
doped SrTiOs ijn^ BM ). fc" and k 1 ' represent the parallel 
and perpendicular plans to the direction of the compressive 
100 strain resulting from Vq. 



Direction 



m CBM 
m he 



DC B M 
m he 



k 1 - 



jfeU 
k 1 - 



r -> x(ioo) 

7.3 
7.3 

T -> M(110) 

1.1 
1.1 

r-»i*(in) 

0.9 



0.5 
11.2 



1.1 
1.1 

0.9 




• • $ • 

• • • • 

FIG. 5. Isosurface of the highest occupied orbitals viewed alonj 
VS (left) VSi (right). 

free electron mass (m^ e /m e ) of the lower band or heavy 
electron (he) band are isotropic and experience a decrease 
from the bulk value of 7.3 to 6.8 in the T — >-X direction, 
bringing the m* he /m e ratio closer to the experimentally 
measured values [3] of 6-7.1 for 15% La doping. A smaller 
decrease is observed along the other high symmetry di- 
rections. 



IV. SUMMARY AND CONCLUSIONS 

The structural and electronic properties of point de- 
fects Vq, Kp i? Vg° r and La^y in SrTiOs have been com- 
puted using the HSE screened hybrid functional. The 
crystal field splittings of the conduction and valence 
bands, resulting from the atomic relaxation around the 
various defects, were evaluated. HSE is known to give an 
accurate description of band gap and VBW of SrTiOs, 
which leads to the assumption it will perform also well 
the defect energetics, as we have shown here. In fact, the 
location of the neutral defect bands in the band gap of 
STO calculated with HSE/SZVP do not suffer from the 
band gap underestimation problems displayed by semilo- 
cal functionals, thus negating the need for further correc- 
tions. The wavefunction is localized around Vq and 
as shown form the isosurface of highest occupied orbitals 
in figure 5 and delocalized for the defect states that al- 
most or completely overlap with host bands like Vg° r and 
Lasy. Our calculated defect bands positions are repre- 
sented schematically in Figure 7; we present the results 
this way for ease of use and comparison to subsequent 
works. The location of the defect level in the band gap 
of STO indicates that Vq is probably at the origin of the 
blue luminescence of STO, and serves as a double shallow 
donor under thermal equilibrium. This defect level dia- 
gram might serve as a guideline in the interpretation of 
photoluminescence experiments. [1, 3] For SrTiOs doped 
with Vq, the calculated conduction band electron effec- 
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the 100 direction showing (a) orbital localization around (a) 




r X M r R X|M R 

K-Path 



FIG. 6. (Color online) Band structure for spin up (Top) and 
down (Bottom) electrons in a 2x2x2 SrTiOs supercell doped 
with La. The Fermi energy Ef is indicated by a solid black 
line. For spin up electrons, one of the triply degenerate bands 
in the CBM is populated (shown in red) while the remaining 
empty ones are shown in blue. 
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Finally, with its two fixed and unmodified parameters, 
the HSE functional gave a reliable description of the elec- 
tronic structure for STO and its defects, in agreement 
with the findings in Ref. 12 and 51. It would be very in- 
teresting to extend the present study for larger supercells 
followed by the a finite size scaling. This would enable as- 
sess the robustness of the present agreement we get with 
experiment and allow a better comparison by performing 
double doping with La and Vq. 
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